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Abstract 

We have considered a semi-infinite crack embedded in a transversely isotropic medium and 
studied two special cases, one, in which the axis of symmetry is normal to the crack face and 
the wave incidence is arbitrary and another, in which the axis lies in the crack plane normal 
to the edge and the incident wave vector is also normal to the edge. The problem is of interest 
in Non-Destructive Evaluation, because austenitic steels that are found in claddings and other 
welds in the nuclear reactors are often modeled as transversely isotropic. In both of cases, we 
have expressed the scattered field in a closed form and computed the corresponding difi^raction 
coefficients. 
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1. Introduction 

The main aim of this article is to evaluate — in two special cases — diffraction coefficients 
of a semi-infinite planar crack embedded into a TI (transversely isotropic) solid. TI is 
an anisotropic medium which is invariant under the rotation around a symmetry axis, 
that is, a TI medium exhibits isotropy in the plane perpendicular to that axis. The 
symmetry of this kind is ubiquitous in nature and manmade materials and TI solids 
have been widely studied e.g. in crystallography, seismology and NDE (Non-Destructive 
Evaluation) . 
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The diffraction coefficients relate the far field amplitudes of the corresponding waves 
diffracted by a straight crack edge to the amplitude of an incident plane wave. The rays, 
along which energy propagates and the wave fronts, which separate the disturbed regions 
in space from undisturbed, are the basic concepts of the so-called ray theory. In the far 
field (high-frequency) approximation, this provides a convenient description of the wave 
phenomena such as propagation, reflection and refraction (see e.g. [2]). 

It is well known that in an isotropic solid, the fronts of waves radiated by a point source 
comprise three concentric spheres. At each moment in time, the sphere of the largest ra- 
dius is the region covered by the disturbance due to P (compressional or longitudinal) 
mode of propagation while the remaining two spheres are coincident and describe the re- 
gion disturbed by S (shear or transverse) waves. The latter degeneracy is visually resolved 
by identifying two orthogonal components of polarization, which are named the SV and 
SH modes, with V standing for vertical and H for horizontal. The rays are normal to 
wave fronts. The nature of wave propagation in anisotropic elastic solids is significantly 
different, since they support three distinct types of elastic waves, qP, qSV and qSH, with 
q standing for quasi, none of which is in general purely compressional or shear. The rays 
are not normal to the wave fronts either [3] 

A well known extension of the ray theory, GTD (the Geometrical Theory of Diffrac- 
tion), suitable for description of diffraction in the far- field approximation was first in- 
troduced by Keller [4]-[6]. He considered diffraction of the plane scalar wave by an edge 
embedded into an isotropic medium whose radius of curvature is much larger than the 
wave length. According to GTD, a plane wave incident on such an edge produces a cone 
of diffracted rays. The apex of the cone coincides with the point of diffraction; it is cen- 
tered on the straight line which is tangent to the edge at this point, and the cone's solid 
angle is determined by the Snell's law. Moreover, according to Keller, GTD applies even 
when the incident wave is not plane, but at the point where it strikes the edge, the radii 
of curvature of the wave front are much larger than the wave length. Special cases which 
can be treated using the high-frequency approximation but for which GTD fails can be 
described with the so-called uniform GTD which lies outside the scope of this article. 

In [1] Keller's GTD was extended to elastodynamics and derived formulae for the 
diffraction coefficients for a semi-infinite planar crack in an isotropic; solid. They achieved 
this by reformulating the original problem as a Wiener-Hopf matrix functional equation 
in a complex Fourier variable. The problem was solved by the Wiener-Hopf factorization 
technique (see e.g. [15]) and explicit formulae for the diffraction coefficients were derived 
using the method of steepest descent. The latter provides an explicit construction of the 
diffraction cones as well as complex amplitudes along the rays. 

For a general anisotropic solid, a semi-analytical approach to the problem was previ- 
ously developed in [13]. The authors have reduced the problem to a Wiener-Hopf func- 
tional equation which, in general, has no known analytical solution and have used a 
numerical scheme to factorize the underlying Wiener-Hopf matrix. In [14], they gave a 
description of numerics and showed the dependence of the magnitude of the backscatter 
qP — qP diffraction coefficient on the observation (incidence) angle. 

Here we aim to produce a procedure for calculating some diffraction coefficients for 
semi-infinite planar cracks in austenitic steels. The problem is of interest in NDE, because 
austenitic steels are found in claddings and other welds in the nuclear reactors. It is well 
known that the austenitic steel can be modeled as a TI material (see e.g. [12]). A simple 
case of normal incidence in a TI material that supports three convex slowness surfaces (see 
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below) has been considered before [7]. We address a more challenging oblique incidence 
case and deal with an extra complication due to the fact that in the austenitic steel, one 
of the slowness surfaces has inflections. 

The article is organized as follows: in Section 2 we state the problem in terms of partial 

differential equations and boundary conditions. In Section 3 we discuss the corresponding 
transfer tensor, that is the free-space Green's tensor in the Fourier domain. In Section 4 
we reformulate the problem in terms of an integral equation (Green's formula). In Section 
5 we reformulate it again as a 3D functional equation and solve this equation analytically 
in Section 6 for the case of the symmetry axis perpendicular to the crack plane and in 
Section 7 for the case of the symmetry axis lying in the crack plane perpendicular to the 
crack edge. In Section 8 we calculate the corresponding diffraction coefficients. 

2. The problem statement 

Let the medium be a homogeneous elastic solid governed by the Hooke's law 

4(x^t)=4.«e«(x^t), (1) 

where the superscript is used to denote the dimensional quantities (whenever the non- 
dimensional versions arc also used), x"^ is an arbitrary point in space, t is time and 
everywhere a'^{x.'^,t) and e{x.'^,t) are respectively, stress and strain tensor of the sec- 
ond order while i,j,k,l = 1,2,3 are indices of the tensor component corresponding to 
any three dimensional Cartesian coordinate system. Here and everywhere below, unless 
otherwise stated, we employ the summation convention over the repeated index. 

As already mentioned above, a TI material has one axis of symmetry. It is well known 
that the corresponding stiffness tensor c"^ involves five unknowns [3]. Also, the stress 
tensor a''(x'^,t) and strain tensor e(x'^,f) are both symmetric, which allows us to re- 
duce their order by using the so-called Voigt notations: Introducing the medium Carte- 
sian coordinate system {81,82,63}, with the 83 axis running along the symmetry axis, 
we replace the tensor cr'^(x'^,t) by the six dimensional vector whose the components 

1 (x'' , t) , a^2 (x'' , t ) , erf 3 (x'^ , i ) and a^^ (x'' , t) , ^ (x'' , t ) , crfa (x'' , t) are called the normal 
and shear stresses, respectively. Then the Hooke's law (1) takes the form 
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where = {Af-^ — Af2)/2. It follows that any TI medium is characterized by five 
independent elastic moduli Af^,Af2, Af^, Af^ and Bf as well as density p. 

To continue, in the absence of body forces, the elastodynamic equation can be written 

as 
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■ a\^,t) = pd^y{^,t), (3) 

where the nabla operator = {df,d2,d^), di is the partial derivative with the re- 
spect to time and u'^(x, t) is the displacement. Let us further assume that the medium 
is irradiated with an incident wave u''^'"^(x, i), which is harmonic, plane and has amph- 
tude C/q . The corresponding stress, strain and displacement fields can be written as 
a<-"'\-}iJ) = crf,P^^(x)cxp (-ioit), g(-"=)(x,f) = eif^(x)cxp {-lojt) and u'^(-"=)(x, t) = 
UTO'"^''(x)exp {—iujt), respectively, where the subscript m denotes functions whose vector 
arguments are expressed in the medium coordinates. Below we simplify the presentation 
by dropping the factor cxp(— icjf) everywhere. A further simplification can be achieved 
by non-dimesionalizing all physical variables, except uj and t which do not feature below, 
using the material density p, the S wave speed along the symmetry axis cq (where a 
degeneracy takes place and there is no distinction between the qSV and qSH modes), 
reference wave number ko = u>/cq and the amplitude of the incident wave Uq, that is by 
introducing non-dimensional variables 

A'^ nd 

pci p4 p4 

k = A:o-V, u = (;7o'')-^u^ 

x = A;ox'', W = kQ^V'^. (4) 



The dimensionless form of the reduced elastodynamic equation is 

V • cr + u = 0, (5) 



where u = am(x) and u = Um(x) . Since for small deformations, the components of 
strain tensor e = Cmlx) are defined by 



\{djUi + diUj) , (6) 



the dimensionless form of the Hooke's law (1) allows us to relate the stress tensor to the 
displacement using the formula 



<Jii=^fk{^)uk, (7) 



where S(V) is a diflferential operator with the elements E^^^, S^^^ and TS^\ which are 
the 3x3 matrices 



4 







^1292 


^13^3 


S(i)(V) = 




Bsdi 







Bids 





Bidi 




Bsd2 







S(2)(V) = 




^11(92 









Bids 


Si 92, 




Bids 





Bidi 


S(3)(V) = 







Bid2 






^1352 


Assds 



(8) 



Note that the stress tensor is symmetric, (7^ = Gji, and therefore, the divergence of 

the stress tensor can be written as 



V • cr = £(V)u, 

where in the medium coordinate system the operator L = i(V) is given by 
iii(V) = Aitdl + Bsdl + Bidl 



(9) 



L22(V) = Bsdi + Audi + Bidi, 

^3' 



J^33(V) = Bidl + Bi^l + Assdl 



j=l,2. 



il2(V)=L2l(V) = (All - 53)91%, 
i,3(V) - L3j(V) = (Ai3 + Bi)djds, 

It follows that elastodynamic equation (5) can be rewritten as 
Lu + u = 0. 



(10) 
(11) 



Let us now assume that the medium contains a semi-infinite planar crack with non- 
contacting faces (see Fig. 1). Then the problem (11) should be supplemented with the 
boundary conditions: Firstly, we assume both crack faces r+ and F" to be traction-free 
(see Fig. 1), so that we have 



r± 



= 0, 



where the traction t = t(x) is a vector defined by 
t = t> ■ a 



(12) 



(13) 



with u -the inner normal to the upper crack face (pointing into the solid). Note that 
using (7), we can write 



t = 5„(V)u, 

so that t is related to u via the displacement-traction transfer operator 
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(14) 



(^™),fe(V) = ;/,EW(V), j,fc= 1,2,3. (15) 

Secondly, at infinity we assume the radiation boundary condition in the form of the lim- 
iting absorbing principle (see e.g. [9]). This implies that once a small positive imaginary 
part is introduced into a wave number, all scattered waves decay at infinity. Thirdly, at 
the tip of the crack, we impose the so-called tip condition, requiring the elastic energy 
of the total field to be bounded (see e.g. [10], [11]). Combined with (11) this reduces to 
the requirement 

t = Oir~l/'), (16) 

where r^^p denotes the dimensionless distance from the point of observation to the tip of 
the crack. 




Fig. 1. The geometry of the problem: the half-plane crack with non-contacting faces r+ and F , the 
medium coordinate system {61,62,63} and the crack coordinate system {61,62,63}. 



3. Transfer Tensor 



Let us introduce the Green's tensor, u'^ — m^(x), the solution of the system of elasto- 
dynamic equations 

L,,(V)uffe(x) + uf,(x) = -<5(x)<5,fc, (17) 

where 6{k) is the delta function and 6ik is the Kronecker delta. Taking the triple Fourier 
transform, which is denoted everywhere by the hat ^, Eq. (17) gives us 

Lr..mf,,{^)-ug{i)=I, (18) 
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where / is the 3x3 identity matrix and as above, the subscript m is used to denote a 
function whose vector argument is expressed in medium coordinates. The solution of (18) 
is called the transfer tensor. 

To continue, the determinant of the so-called Kelvin-Cristoffel matrix Lm{^) — ^miOI 
can be expressed in the form 



where we employ the notations 



(19) 



Ai_2(|) = [And + BiCl - A™(0][i?id + ^33^1 - A„(^)] - (A13 + i^O^eiCa', 
A3(0 = S3el+BiC|-A™(0, (20) 

with d = $1 + ^2- Thus, eigenvalues of Lm{^), that is zeros of Ai_2 and A3 are 

aL^^(0 = 

J[(Si+An)Ci + (Bi+A33)C3- 



(-1)'^^ { [(An - Bi)?i + (Si - Ass)es? + 4(Si + Ai3)^CiCl} 
and the corresponding unit eigenvectors are 

where no summation is implied and we use 



/3 = 1,2, 



(21) 



(22) 



(Bl + Ai3)66 

(Bi + Ai3)66 



Note that if ^2 0, as ^1 0, we have 



P^H0-e3, P^^(0-ei, p^^(4)-e2. 



V(2) 



(3), 
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(23) 



(24) 



The unit eigenvectors of L„i{^) form an orthonormal basis, which can be used to expand 
any function, that is, we can write 



(25) 



where the outer product ab is a tensor with components (ab),fe = aibk- Here and every- 
where below, unless stated otherwise, we imply summation over the repeated superscript 
/3=1, 2, 3. 

Substituting (25) into (18) and using the definition of eigenvectors, we find that the 
vectors a*^'^) satisfy the matrix equation 



[A(f)(0-l]a('^^plf(0-/- (26) 
Dot-multiplying (26) by Pm\ (3 — 1, 2, 3 from the right, we obtain the transfer tensor 



(27) 



3.1. Slowness surfaces and wavefronts 

In view of the above, the TI media are conveniently described by the so-called slowness 
surfaces generated by the slowness vectors [cm\n)]^^n , where c is the wave speed and 
the unit vector n indicates the direction of wave propagation [3]. The slowness curves, 
cross-sections of these surfaces with ^2=0 are shown in Fig. 2 (a). The slowness surfaces 
can be used to construct the wavefronts of waves radiated by a point source at the origin. 
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Fig. 2. (a) The slowness curves and (b) wave curves for ^2 = 0. When the wave curve is rotated around 
the X3 axis, points A circumscribe the so-called cuspidal edges. Points B are called the conical points. 
Key: qP (solid line), qSV (dashed line) and qSH (dotted line). 

Indeed, the front of each infinitesimal plane wave of type /3 radiated by this source in 
the direction of a unit vector n moves to, after a unit of time, the location x such that 
we have 



n.x = cL« 



(n). 



(28) 



The envelope of the above plane fronts is obtained by differentiating the wavefront equa- 
tion (28) with respect to n. Therefore, the front of the wave radiated by the point source 
is generated by the rays 



dn 



(29) 



The wave curves, the cross-sections of the above wave fronts with ^2 = are shown in 
see Fig. 2 (b). Note that the rays indicate the direction of the energy propagation and 
are perpendicular to the slowness surfaces and not the wavefronts (see Fig. 3 (a)). 
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4. The reciprocity theorem (Green's formula) 



Let us introduce the Cartesian system {ei, 62, 63} associated with the crack, such that 
its 63 axis is perpendicular to the crack plane, 62 runs along the crack edge and ei is 
the inner normal to the edge. Note that there exists a matrix of coordinate rotation Q, 
such that = Q^, where the superscript denotes a transpose, and the medium 
coordinates Vi of any vector v = k, n, x, ^ etc. are related to its crack coordinates Vj via 



1,2,3. 



(30) 



Below, the bold symbol v is used to denote the triplet of coordinates of vector v in the 
primed (crack) coordinate system. When working in the unprimed (medium) coordinate 
system the symbol v is retained to denote the triplet of the corresponding coordinates. 
Since the medium is isotropic in the eie2-plane, we can always rotate the medium coor- 
dinate system around the axis of symmetry e3 so that 63, the normal to the crack plane, 
lies in the e2e3-plane and therefore, ei lies in the crack plane 6162. Then without loss of 
generality, the matrix of coordinate rotation is 



Q 



cos (^0 sin (^ocos ^0 sin (^osin 60 
—sin (^0 cos ^ocos 6q cos (^osin 6q 
—sin ^0 cos 6q 



(31) 



where 9q is the angle between the symmetry axis and normal to the crack, and 0o is the 
angle between ei and ei. 

We can now use the reciprocity theorem (see e.g [1]) to reformulate the problem stated 
in Section 2 as an integral equation for the scattered field (otherwise known as the Green's 
formula). We start by decomposing the total fields u and a as 



u = u'" + u'" 



(32) 
(33) 



with u™ and cr™ representing the respective incident displacement and stress fields. The 
incident displacement field can be written as 



u*"(x) =p(^)(k*")e*'"-'^, 
where we use the notation 

u-(x) = uj;:(x) ■ 

x=QJ X 

the polarization vector p^'^^ = p^^-'(4) is 



{0)1 



(34) 
(35) 

(36) 



the wave vector is k'" = fc^'^^h'", where n*" is the propagation direction of the incident 
wave, given in the primed (crack) coordinates, and the wave number fc(^) = a;/c(^)(n'") 
is 
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A;'" = [AW(n'")]-V2^ (37) 
where similarly to (35) and (36), we use the notation 

A(^HC)=A(^^'(€)|^^^^^. (38) 

We note that the incident stress field gives rise to the incident tractions 

t'"(x) =t('^)(k*")e*'"-^, (39) 
and the boundary condition (12) can be rewritten as 

t*^(x) =-t'"(x) . (40) 



r± 



r± 



The amplitude of the incident tractions t('^)(k™) is related to the amplitude of the 
incident displacement field p^''^ via 

t(«(0 = i5(^)p(^)(0- (41) 

The latter relationship is obtained by substituting (34) into (14), applying the triple 
Fourier Transform to the result, changing from ^ to ^ and finally using (15) and the 
notation 

Si^) = Sm{$)\^^^^^. (42) 

As before, the presence (absence) of subscript m means that the function has the vector 

argument expressed in the medium (crack) coordinates. 

By taking into account that both crack faces are traction-free, we then write the 
reciprocity theorem for the scattered field u*'^ as 



u^^(x) = - / dy^ dm t^(i - y)Au^^(y), 

Jo J-oo 



(43) 



where t'^ = u-a^ is the Green's traction tensor, u = 63, cr^^, = Ej^^(V)u^ is the Green's 
stress tensor and Au*^(fi,f2) is the so-called COD (Crack Opening Displacement) (see 
e.g [1]), which is defined as 

Au^"(f 1, f 2) = u^"(fi, f 2, 0+) - u^=(f 1, f 2, 0-). (44) 

Above, 0+ refers to the upper face of the crack and 0— to the lower. 

Let us note that the scattered field is invariant with respect to translations along 
the edge of the crack. Therefore, all the fields u, t, a and Au**'' have a common factor 
exp(— 1^2^2)1 where we have 

6 = -kT, (45) 
and all can be factorized as follows 

v(x)=V(fi,f3;<e2)e-'«=*^ (46) 

Here and everywhere the argument ^2 is separated from other arguments by a semi- 
column to emphasize the fact that in the problem under consideration it is just a fixed 
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parameter. Also, the bold capitals are used everywhere to denote the preexponential 
factors of the quantities denoted by the corresponding lower case letters. Then the integral 
equation (43) can be rewritten as 

u^^(ii,f3;6) = 

- dm dy2t'^{ii-yuX2-y2,X3)AV{yi;i2)e'^'^^'-^'\ (47) 

Jo J-oo 

where the components of vector functions U*'^, AU*'^ and tensor t'^ are all given in 
medium coordinates while the crack coordinates are used to represent the components 
of vector arguments. This choice leads to simpler formulas. 



5. The 3D functional equation 



Noting that the t/2— integral in (47) is actually a single Fourier transform, which is de- 
noted everywhere by the bar ", and applying the convolution theorem to the yi— integral, 
the single Fourier transform of (47) in xi gives us the vector functional equation 

U''(a,i3;6) = -F(^'i,i3;6)AU''(6;6)- (48) 

Above and everywhere below, the tilde ~ denotes the double Fourier transform The 
tensor t'^{^i,xs;^2) is the inverse Fourier transform in ^3 of the traction tensor ?'(^), 
which is given by 

1-A('3)(6 ^ ' 

This means that we can write 
~ 1 f°° 

t''(6,^3;6) = ^y_ i^(0e-^«3^3rf6- (50) 

Note that as ^3 00, t^{t) = Therefore, above we can apply the Jordan Lemma 

if for fa > 0+ the contour of integration is closed in the lower ^3-plane and for X3 < 0— , 
in the upper half-plane. In both cases considered below, for each ^1, inside the chosen 
contour, the Transfer Tensor has two poles ^"■'^(C'l; C'2) a = 1,2, which are roots of the 
quartic expression in (20) with the corresponding A*^")(^) = 1, and one pole 'f3^'^(Ci; C2), 
which is a root of the quadratic expression in (20) with = 1. Above, the top sign IS 

chosen when X3 > 0+ and bottom when £3 < 0—; and for the cases under consideration, 
Cl"^'^(a;6) = TC3"''(6;6); where the functions 6) are defined in (59) and (99), 

respectively. The real parts of functions (59) are shown in Fig. 3 (b). When $l"^(^'i;6) 
run along slowness surfaces the functions have no imaginary parts. Fig. 3 (b) shows that 
a = 1 always describes the qSV mode, a = 3 -the qSH mode and a = 2 — the qP or 
qSV mode, depending on ^i. The functions (99) exhibit analogous behavior. 

On adding a small negative imaginary part to all c^"^ (a small positive part to the 
wave numbers), each real pole moves away from the real axis in the a— plane, with the 
sign of the resulting — Im ^3"^^(a, c^"^) the same as the sign of the partial derivative 
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1 

Fig. 3. (a) The qSV energy flow (the mean Poynting) vectors P = — ilm [ojcr'' • (u**)*], where the star 
denotes the complex conjugate, (b) The real parts of three roots ^3^''(^i;^2) (dashed line), $3^''($i;^2) 
(solid line) and ^'3^^ (^'i ; ^2 ) (dotted line) imposed on the slowness curves in Fig. 2 (a). When lying on 
the slowness curves the roots are real. When sliding of these curves at the branch points ^1 = ki,K2,kz 

or K4 the roots arc complex. 

5^3"^^(^i, c^"^)/5c^"). In other words, the resulting — Im ^3"^^(Ci) c^"^) has the same sign 
as the vertical component of the gradient to the slowness surface at ^1 (see Figs. 2 (a) 
and 3 (b)), that is the sign of the vertical component of the energy flux (see Fig. 3 (a)). 
Since P (4) in (49) is a rational function of ^3 we can now evaluate the integral (50) using 
the Cauchy Residue Theorem to obtain 



«'3=4"'*(€'i;«'2); A(")(^)=i 



where the subscript ^3 indicates differentiation with respect to ^3 and the resulting plane 
waves satisfy the radiation condition at infinity. For the two cases considered in this 
paper the denominators in (51) are specified below in (124) and (125), respectively. We 
also note that i'^(^i, 0— ; ^2) — ^'^(Cii 0+; ^2) = I, the identity matrix. This can be best 
seen by using (48) to evaluate U*'=(^i, 0+; ^2) - U^'=(6, 0-; 6)- 

In order to utilize the boundary condition (12) in a straightforward manner, we apply 
the operator — iS'(^i, ^2, 1^3) to Eq. (48) and then write the resulting equation for the 
crack face X3 = to obtain 

T''(6,0;6) = -T^(a,0;6)AU''(6;6), (52) 

where 5(^1 1^2,1^3) = S{!^), with 18/8x3 substituted for ^'3, and we have 



T (4i,a;3;^2) = ±1 



A^^(^) 



. (53) 

«3=4"'^«i;€2); A(»)(^)=l 



Note that on the crack plane, T*'^ must be a continuous function of x, so that T (^1 , 0+; ^2) 

T (^1,0—; ^2) and therefore. t'~^{$,i,0;^2) is well-defined. This can be verified indepen- 
dently by expanding T^{^), the inverse Fourier transform in ^3 of (53), into the Laurents 
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series, to verify that for any crack orientation, its components exhibit the following be- 
havior at infinity: 

7^(4) = constant matrix + 0(4-), as ^'3 — > oo. (54) 
Then t^{^i,0;^2) is well-defined via the following consideration: 

T^(6, 0+; 6) - 0-; 6) = lim / ?^(^) rf^s = 0, (55) 

where the integration contour is a circle of radius R, circumscribed anticlockwise and 
centered at the origin of coordinates. In view of the above, below we use the simplified 
notations 

T-(fi;6) = T-(fi,0;6), t«(6;6) = r°(6,0;6), (56) 

We can now make use of the boundary conditions (40) and apply the Fourier Transform 
to (39) to put the functional equation (52) in the Wiener-Hopf form, 

^ 6 +T?iO -^^°"(^'^^^'^) -^"(^'^'^'^)^^^(^'^'^'^)' (57) 

where AU'*'^(^'i; ^'2) and T*'^ (C'ljC^) are both unknown. Since AU'*°(f 1; ^2) = and 
T'^"^^ {xi; X2) = for .x 1 < and .x 1 > 0, respectively, their Fourier transforms AU (^1; ^2) 
and T"^ (^i; ^2) are analytic in the upper and lower half of the complex ^i-plane, respec- 
tively. This can be seen by studying the corresponding single Fourier integral. Above, we 
also use the fact that we have 



L 



gi[€i+fer]c6i^^^ = , . (58) 

a + M" + «o 

6. Solution of the 3D functional equation when the symmetry axis is 
perpendicular to the crack plane 

6.1. The poles and branch points oft'^{^) 

When the symmetry axis is perpendicular to the crack plane, 4>o = 0o = 0, Q = I and 
thus, all the ^-components coincide with the corresponding ^-components. As mentioned 
above, the functions ±^'3"^(^'i; ^'2), a = 1, 2 and ±C3^^(^'i; C2) which are the poles of the 
tensor are the respective roots of the quartic and quadratic expression in (20), with 

A(")(^) = 1, so that we have 

^?^(ei;6) = 61(51; 6) +62(6; 6), 
5?^ (a ; 6 ) = isi (6 ; 6 ) - 62 (a ; 6 ) , 

(6; 6) = (i?r'i?3)^/'74(6; 6), (59) 

with 
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^3n(^i;^2) — ^ 



2(SiA33)l/2 

[Al^ + 2BiAi3 - AnA33]il +Bi + A33 - (-l)"2Bi(AiiA33)i/27i72 1 



76-3n 



1,2, (60) 



{}^/^ denoting the principal value, 



7o=7o(Ci;6)-(«:^ + 4')'/' 

le = le{iiU2) = [{ne + i0f-i!]'^^ £ = 1,2,3,4, (61) 
with the branch points ii/to and ±Ke, £ = 1,2, 3, 4, where we have 

K4 = {B3'-ilY'^ (62) 

and 

Ci = £)-i{2(A33-Bi)i/2(s^ + ^13)1^23 - ^11^33 + 5i(An + 2Ai3 + ^33)]^/' - 
(-l)^[Ai3(2Bi + Ai3)(Si + ^33) + ^33(2Bf + AnSi - A11A33)]}'/', £ = 0, 3, 
D = {(^11^33 - A?3)[(2Bi + A^sf - ^11^33]}^/'. (63) 

It can be checked, e.g. using MATHEMATICA, that when n = 2 (n = 1), ±iKo (i'ts) 
is a zero of the corresponding numerator in the curly brackets in (60). Therefore, the 
square roots of the expressions in the curly brackets never vanish. Note that all branch 
points are either real or purely imaginary and for the austenitic steels that are described 
e.g. in [8] and are considered in this article, direct evaluation shows that we have the 
following inequalities: 

Re Ki < Re K2 < Re /ts < Re K4, 

Im K4 < Im K3 < Im K2 < Im ki < kq. (64) 

Therefore, the phases of the principal values of ^32 and ^31 are determined by the position 
of ^1 on the complex ^i-plane relative to K2 and K3. For a fixed values of ^2, the real parts 
of functions ^s"'', a = 1, 2, 3 described by equations (59) are represented in Fig. 3 (b). 



6.2. Decoupling the normal COD from the tangential COD 

In the case under consideration, t'^ = cr'.^ and the crack and medium coordinates 
coincide. Therefore, in view of (42) the transfer operator S{^) is 
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Bi6 BiCi 

S{i)= Blis 

^136 ^336 

It can be checked, e.g. using MATHEMATICA, that we can write 

M3(a;6) 



(65) 



T^(6;6) 

where we have 

M3 



^ / -433^/2 72(6; 6) .. , X 



71(6; 6) 



and 



with 



^'(6;6) 



6 6 

6 -ii 



(66) 



(67) 
(68) 

(69) 



M(6;6) = 



Ml (6; 6) 
M2(a;6) 



(70) 



Ml (6; 6) 



ii?(6;e2) 



4^3361(^1:6) 
M2(6;6) = ^(-SiB3)'^%4(6;6), 

and the Rayleigh function is defined by 

R{ii;i2) = (A?3 - AnAssM + il) + ^33 + (An^33)^/^ ^^Sf'ti - 

72(6; 6) 

Then substituting (66) into the 3D vector functional equation (57), the latter decouples 

into the scalar equation describing the COD component normal to the crack face and 2D 
vector functional equation describing the tangential COD components. 



(71) 



(72) 



M3(6;6)AC/r(6;6) = i , - Tr(6;6) 



(73) 



and 
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Mfe;e2)AZ^^^(6;6) = i . ^ S J -T^^'(6;6), (74) 

where the calligraphic script is used to denote two dimensional vectors expanded in the 
eigenvectors of r'"'^'', so that for any vector A. = Atl \ T^'^ or f^^'^ we have 

Miv,i2) = 6), i, k = i,2, (75) 

where we remind the reader that the repeated index summation is applied to fc = 1, 2. 

6.3. Solving the scalar functional equation for the normal COD 

First, let us consider the scalar equation (73). It can be solved using the Wiener-Hopf 
technique by introducing the following factorization of the tangential traction ii^: 

M3(6; 6) = (6; 6)m3 (a; 6), (76) 

where +(— ) means that the function is analytic in the upper (lower) half of the complex 
^1 -plane. In order to effect this factorization, we first note that there exist values = 
±Ki{, kr = {k\ — Ci)^^^ such that 

i?(±Kfl;6)=0, (77) 

which arc called the Raylcigh poles. Above, kn = cq/cji, where cr is the speed of the 
surface Rayleigh wave. Only the pole —kr corresponds to the outgoing wave. Then 
the considerations and notations introduced in Appendices A and B allow us to define 

Mf(£.;{2), 

MfK.;4) - (j^f Jf'"!:!^ ., ('''-^^-). P8) 



and therefore, A*3'(^'i;^2), 

where throughout we use the notations 

7o'=(ei;e2) = (iA^o±a)'/', 

7^(e'i;6) = («^ ±6)'/', ^= i,...,4. (80) 

To continue, as mentioned at the end of Section 5, AU''^(^'i; ^'2) is analytic in the upper 
half of the ^i-plane. Hence Eq. (73) can be rewritten as 



-/x^(6;6)A[7r(6;e2) + i- — ^'^^^"^^ , =tT''~, (8i) 

(ei + A:r + ^o)M3-(-fcr;6) 



where the superscript """^ stands for modified and we have 

16 



=TO0d- 

3 



Ma" (a ; 6) M3 (-fcf; 6) Ma" (a ; a) a + M" + ^0 



(82) 



Note that the modification leading to T^'"^ has been introduced to assure that the 

left-hand (right-hand) side of (81) is analytic in the upper (lower) half of the complex 
^1— plane. This can be true only if both sides are one and the same entire function. 
Furthermore it is easy to check that at infinity we have 



nt{ii-A2) = o{i\'\ 



(83) 



1 /2 

Also, the tip condition (16) implies that AU'*'^(f i, f 2) = 0(fj ), which in its turn 
implies 



(84) 



It follows that as ^1 00, the left-hand side of Eq. (81) and therefore, its right-hand side 
are both functions of order 0([^i]~^), obviously bounded. According to the Liouville's 
Theorem, any such entire function is in fact zero. Therefore, Eq. (81) implies 



Af/a (a;a) 



t^3''^(k^") 



(a + M" + ^o)m^ (a; a)M3 a) 

6.4. Solving the 2D vector functional equation for the tangential COD components 



(85) 



Let us now turn to the 2D vector functional equation (74). The matrix M can be 
factorized. 



M(a;a) = M+(a;a)M-(a;a), 

so that the matrix 

M^(a;a) 

is analytic in the upper half of the complex ^1— plane, and matrix 



(86) 



(87) 



M-(6;a) = 



Mr (a; a) o 
M2(a;a) 



(88) 



is analytic and has no zeros in the lower half of this plane. The factorization of ^ii{ti')£,2) 
has been described in the previous section, and /i2 can be readily factorized into the 
factors 



M^ 



(a;a) = (-^)''V(a;a). 



(89) 



Below we also use the matrices which are respective inverses of . Now, multiplying 
(74) by matrix N~{^i;£^2) gives us 
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M+(6;6)AZ^''(6;6) = iiV-(6;6) 



It can be transformed into 



(90) 



where we have 

■r=mod— , ' ' 

^ (a; 6) = 



6 + ^i" + iO 



(6; 6), (91) 



A^-(a;6)T^ (6;C2) + i[A^-(-fcr;6)-A^-(ei;6)]^^^r^, 



iO 



(92) 



so that the left-hand (right-hand) side of (91) is analytic in the upper (lower) half of 
the complex ^i— plane. This means that both are one and the same entire function. It 
is easy to sec that as oo, this function is 0(1) and therefore, bounded. According 

to Liouville's Theorem, any such entire function is in fact a constant. Let us call it g. 
Then multiplying both sides of (91) first by (M+)-i(^'i; ^2) = -/V+(Ci;6) and then by 
P-Hiui2) = P{iui2)/{ei+^i) we obtain 



/ -iA^-(-fcr;6)P(C'i;6)f'^Hk") 

TT~{ii-A2)i = \p-\ii;i2)M-{irA2) ■ 



+g 



(- 



[N-{-kr-A2) - iV-(a;6)]P(C'i;6)f(''nk*") 



+ g 



i = l,2, 
(93) 



where wc have employed the formula (41). In order to find g, we note that in view of 
(69); the plus function AU*'^(^i ; ^2) has an apparent pole at ^1 i^2; similarly, the 
minus function T (^'i;^^) has an apparent pole at —1^2- Since this is impossible, both 
corresponding residues must vanish. This gives us the following linear system for g: 

Pfe;6)^ni6;6)( -'^'<-^yf ^t""''"'"' +«) °°- 

\ -I- -I- ill / 



i6 + fcf + iO 

P(-i6;6)M-(-i6;6)- 

/ -i[A^-(-fcr;6) - iV-(-i6;6)]P(-i6;6)r<^^(k' 
^ -16 + fcf + ^0 



= 0. 



(94) 



Since |P(iC2;C2)| = \P{—''^^2]£,2)\ = 0, only two of the four equations above are linearly 
independent and thus, g is well defined. 
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7. Solution of the 3D functional equation when the symmetry axis lies in 
the crack plane perpendicularly to the crack edge 



7.1. The poles and branch points ofi'^{^) 



When the symmetry axis lies in the crack plane perpendicularly to the crack edge we 
have (po — So — tt/2, 



g = 



1 
-10 
0-10 



(95) 



and therefore, 

6 = -6) 6 



-6. 6 = Ci- 



(96) 



The case of a normal incidence, ^2 = 0, lends itself to an easy analytical treatment. 
Indeed, the poles of i'^{^i;0), are respectively, the roots of the quartic equation 



(SiCf + Ants - l)(^33^f + B,t, - 1) - (Ai3 + Bi)^Cf^3 = 
and the quadratic equation 

0, 



(97) 



(98) 



which correspond, respectively, to the first and second lines in (20), with the correspond- 
ing A'"^ (^) = 1, both rewritten in terms of the crack coordinates. Moreover, Eq. (97) may 
be obtained from the first equation in (20) simply by putting ^2 to zero and allowing An 
and ^33 to exchange places. It follows that we can follow the form of solution presented 
in Section 6.1 and write 



e?^(a;0) = e3i(a;0) + e32(6;0), 



(6:0) = 61(6; 0)- 62(6; 0), 

if\ii;Q) = {B^'B,Y'^^^{i,-Q), 



(99) 



with 



76-3rt 



6n(a;0)- 2(5^^^^)l/2 " 

f [^?3 + 2B1A13 - AnAsi]il + Bi + An - (-l)™2Si (^11^33)^/^7172 l ^2 



76-3n 



and 



7o = 7o(6;0) = (^^ + e?)'/', 



= 1,2,3, 



n=l,2, (100) 



(101) 
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with the branch points ±iKo and ±Ke, £ = 1,2, 3, where we have 



"-^0; '^l ^ ^33 ' '^2 — ^1 J 



— 



(102) 



and 



Ce = 



£>-i{2(Ani3i)i/2(ij^ + A,3)[Als - A^A^s + BMn + 2^13 + ^33)]^/' - 
(-l)^[Ai3(2Bi + Ai3)(Si + An) + Aii(2B? + A33B1 - A^^A^^W^, 

£ = 0,3, 

D = {(^11^33 - Al^)[{2Br + A13)' - A11A33]}'/'. (103) 
The factorization described in Appendices A and B is appHcable provided we have 



Kl < «2 < K3- 

This condition is satisfied by the austenitic steel under study. 



(104) 



7.2. Decoupling the vector functional equations into the scalar functional equations 

In this case, = cr^g ^'^^ the components of the transfer operator ^^(V) are given by 

(^™)ifc(V) =S(^)(V). (105) 
Therefore, in view of (42), the transfer operator S{^) at ^2 = is 



^(a,o,6) = 



-B3C3 
-A11C3 AiaC'i 
B^ii -Biis 



(106) 



It is easy to check that if ^2 = and An and A33 exchange places the Fourier transform of 
the elastodynamic equation is the same as in the previous case. This means that r'~'(^i ; 0) 
has the following diagonal form 



/ii(6;0) 
T^(a;0)= /X2(ei;0) 

M3(a;0) 

with the eigenvalues 

Mi(a;0)-i(-i3iB3)'/'72(a;0), 
M2(6;0) = (^l ^, M3(6;0), 

^^33^ 7i(^i;0) 



(107) 



(108) 



M3(6;0) 



ifi(ei;0) 

4An6i(a;0)' 
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and the Rayleigh function 

R{ii;0) = {Ais - AnAss)i + An + (An^ss)'/'^^^^. (109) 

72(6; 0) 

Substituting (107) into the 3D vector functional equation (57), the latter decouples into 
three scalar equations, 

M.(6;0)Ar7r(6;0) = i T^^^;^^'!^ -T^r (ei;0), z = 1,2,3. (110) 

The equations have the same structure as (73) and therefore, can all be solved using 
the Wiener-Hopf technique by introducing the factorization 

/i,(a;0) = /i+(a;0)Air(ei;0), z = 1,2,3. (Ill) 

On allowing An and A33 to exchange places, the tangential traction /i3(^i; 0) is the same 
as the tangential traction /xi(^i; 0) in Section 6. Therefore, it can be factored in the same 
manner, with Iq, h, Kq and Ki the same as in the Appendices A and B but with An 
standing in place of ^33 and vice versa. 

It follows that the considerations and notations introduced in Appendices A and B 
allow us to define /i^(^i; 0) as 



Ml 



(6;0) = (-^)^%r(6;0) 



33^ 7f(a;0) 



^3^^!' ^ \ Ah An J 7,±(£i;0)i^±ffi;0)^ ^ ' 



73^(6; o)i^f (6:0)' 

where, similarly to (80), we have 



7^(6; 0) = {Ki ± £ = 1, 2, 3. (113) 

It is easy to see from their definition that the + (-) functions above have no zeros in the 
upper (lower) half of the complex ^1— plane. It follows that we have 

AC/r(a;0) = i^ ^ '-^-^ ^ , i = l,2,3. (114) 

(6 + fcr + ^o)M+(ei;0)M-(-fcr;0) 

8. Diffraction coefficients 

Let us consider the equation 

U^^(ii,i3;6) = -^y t''(^i,i3;e2)AU''=(a;6)e-'«^*M6, (115) 

which follows from Eq. (48) via the inverse Fourier transform in ^1. Substituting (51) 
into (115) leads us to 
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u^'=(ii,i3;6) = ±^ 



r 

J — c 



p(o)(^)t(")(^) 



«3=4 '*(«i;«2); A(-)(^)=i ^ 

(116) 



where summation over a is implied and the top (bottom) sign is chosen when f 3 > 0+ 
(^3 < 0-). 

Let us now introduce the polar coordinates f±_ and 9 such that we have 
xi = r_L cos 6, 

X3 = r±sin 9, ^'e[0,27r]. (117) 

When values of k^'^r± arc large, the main contributions to (116) come from the stationary 
phase points and other critical points that are described e.g. in [8]. The method of the 
uniform stationary phase (see e.g. [17], [16]) allows us to treat the situations when the 
critical points coalesce, that is when an observation point lies in a transition zone between 
geometrical regions, but here we treat only those observation points for which the critical 
points are isolated, that is the observation points that lie in geometrical zones. Thus, the 
formulas given below arc not applicable near the shadow boundaries, cuspidal edges 
and conical points of the qSV wave surface or points of tangential contact of the qSH 
and qSV wave surfaces (see e.g. [8] and [20]) The contribution of each isolated phase 
stationary point can be evaluated using the standard stationary phase formula, 



27ri 



1/2 



g{t)e"mdt^g{c) ^ e-/(^), (118) 

where the dot denotes the derivative with respect to the argument (see e.g. [18]); and c 
is such that /(c) = 0. In both symmetric cases considered above, the phase function in 
(116) can be written as 

-6^1 6)1^31 - -r±[eiCos ^-Cl"^(6;6)|sin 9]]. (119) 

Let us call each solution ^1 of the equation 

-cos 9 + (6; 6)|sin 9\ = (120) 

the stationary phase point ^["K Then applying (118) to (116) and multiplying both sides 

of the resulting formula by cxp {ik'!^X2) the main contributions to the scattered field in 
the geometrical zones are given by the stationary points £!f\ so that for each fixed a 
and /3, we have the GTD approximation 

u(")*^/(x) - J-D("^)(k*//)e-*"''-^, (121) 

where using (45), the dimensionless diffracted wave vector is 

i<uff ^ (^■(«),_fc-,_sgn(sin 9)it\i[^^;-kT)); (122) 
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and the vector diffraction coefficients D^"^^ are such that we can write 

D(a/3)(k*//) = _sgn(sin 9) ^- j- • 

27ri|sin^'|^^"^(fcf-^^;-fc™)l 



p(")(k*//)t(")(k*-''^') 



with no summation over a. We remind the reader that the first superscript a = 1, 2 or 
3 describes a diffracted wave and the second superscript /3 = 1, 2 or 3, an incident wave. 
In the case of the symmetry axis perpendicular to the crack plane we have 

A(")f£) = -2i (-4?3 + 2^1^13 - A.MCl - + {B, + A,s)X^-\i) 

«^ (Bi+^n)Ci + (Si+A33)e|-2A(«)(0 

a=l,2, 

xfii) = 2B,is, (124) 

where when |4| = |^3"^(^'i; -^2")!, >^^°'K^) = 1- In the case of the symmetry axis lying 
in the crack plane perpendicularly to the crack edge the above formulas apply for ^2=0 
if All and A33 exchange places, so that we have 



Ai"^(^) = -2e 



{Al^ + 2BiAi3 - ^11^33)^'? - 25^1142 + (Bi + Aii)A(")(^) 



"^^ (Si + ^33)^1 + (Bi + Aii)^'f - 2A(«)(4) 

Q! = l,2 



xf{0 = 2Bsis. (125) 

Let us now introduce the spherical polar angles ip and 1? associated with the crack 
coordinate system. Then each incident unit wave vector can be expressed in terms of its 
medium coordinates as 

n*" = (sin <^*"cos 1?*", cos lyJ*", sin (^*"sin 1?*") . (126) 

When the symmetry axis lies in the crack plane and is perpendicular to the crack edge we 
only consider the incident vectors that are perpendicular to the crack edge, i.e. we assume 
_ qqO Using the transpose of (95), this yields the following medium coordinates of 
the incident unit wave vector n'": 

n™ = (0, -sin cos #") . (127) 

The incident wave number k^^^ and parameter fc" can be found by using the formulae 
(37) and (45), respectively. In each case, the polarization of the incident wave p(^)(k*") 
is given by (22) and (36). 

Thus, each incident wave vector of type /3 = 1, 2 or 3 produces three families of 
diffracted wave vectors of type a = 1, 2 or 3, each covering a portion of the cone surface, 
qP, qSV or qSH, respectively, with each cone having a cross-section in the shape of the 
respective slowness surface (see Fig. 3 b). Each wave vector can be expressed as 
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k*// = -(^(")(^?; -fcf ) cos k^, ^("H^; -^^2") sin 1?), 1? € [0, 27r], (128) 
where in the case of the symmetry axis perpendicular to the crack plane we have 



2Di 



2Di 
1 - Bsik^f 



with 



B3 cos + Bi sin 



Di = B1A33 sin''?? - B4 cos^i? sin^-!? + AnBi cos'^t?, 
D2 = (2An[fcr]'5i - Si - ^11) cos^?^ - (B^k'^f + A33 + Si) sin^^?, 
1)3 = {An[kT? - l){Bi[kT]^ - 1), B4 = Af3 + 2A13S1 - ^11^3 



1/2 



(129) 



i33, 



(130) 



and in the case of the symmetry axis lying in the crack plane perpendicularly to the 
crack edge we can write 



^W(a?;0) 



-D2 + {Dj - ^DiY'^ 
-i?2-(£'i-4£)i)V2 



2Di 

1 



1/2 



1/2 



-,1/2 



Bi COSH + B3 sin^4 



(131) 



with 



Di = BiAn sin^i? - B4 cos^^ sin^ + A33B1 cos*t?, 
D2 = -{Bi + A33) - {All + Bi) sm^-d. 



(132) 



Using the fact that the wave front is a polar reciprocal of the slowness surface, it can be 
shown that in all cases, the polar angles and ^ of the ray and wave vector, respectively, 
are related by the following formula 



,?-tan-fi^). 



(133) 



Similarly to [19], in Appendices C and D we plot for various incident polar angles 
and 'i?™, the magnitudes of the diffracted coefficients D"''^ , where 7, (5 = gP, qSV or qSH 
versus the diffracted polar angle ^ € [0,27r], thus tracing their variation over the corre- 
sponding slowness curves. As expected, all the graphs for the magnitudes D^'''^^''^^^^ are 
similar to the respective graphs for the magnitudes £)SH,sh giy^j^ ^^9]. Two represen- 
tative graphs of the magnitudes of the diffracted coefficients £)('J^''J^) versus the physical 
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e, ° e, ° 

Fig. 4. The magnitude of the diffraction coefficient d'^i^^i^^ versus the physical polar angle d for the 
symmetry axis perpendicular to the crack plane (solid line) and the symmetry axis lying in the crack plane 
perpendicularly to the crack edge (dashed line). The incident angles are 0'" = 90° and (a) = 120°, 
(b) e*" = 150°. 

polar angle G [0, 27r], that is tracing their variation over the corresponding wave curves 
are given in Fig. 4. 

Finally, we note that the GTD approximation (121) breaks down in the so-called 

transition regions between various geometrical zones when the denominator in (123) is 
small or vanishes. The regions center on the light-shadow boundaries, cuspidal edges and 

conical points where, respectively, U^^^ + fc", (ji^^^ \ —^2") ^"^^ A^"'(k*-^^ ) turn to 
zero. When the singularities are few and far between, like in the isotropic case considered 
in [19] or in the cases presented in our Appendices C and D, there is no practical need 
for developing the Uniform GTD that provides a good approximation of the diffracted 
field in transition zones as well as the geometrical ones. However, when cuspidal edges 
and conical points are present the formula (123) ceases to be of much practical use and 
the Uniform GTD description should be used instead of (121) (see e.g. [8]). We illustrate 
the situation in Appendix E but do not plot the corresponding " diffraction coefficients" 
in Appendices C and D. 

9. Conclusions 

We have considered a semi-infinite cr;u;k embedded in a transversely isotropic medium 
and studied two special cases, one, in which the axis of symmetry is normal to the crack 
face and the wave incidence is arbitrary and another, in which the axis lies in the crack 
plane normal to the edge and the incident wave vector is also normal to the edge. The 
problem is of interest in NDE, because austenitic steels that are found in claddings and 
other welds in the nuclear reactors can often be modeled as transversely isotropic. In 
both cases, we have expressed the scattered field in a closed form and computed the 
corresponding diffraction coefficients. 

A simple case of normal incidence in a TI material that supports three convex slow- 
ness surfaces has been considered before [7]. This article and its published version [21] 
addresses an extra complication which arises when a slowness surface contains inflections 
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and in the first of the above cases, the incidence is allowed to be oblique. In future, we plan 
to cross- validate our code with other numerical codes, such as the one reported in [13], 
and also validate it against experimental data — whenever the latter become available. 
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Appendix A. Factorization of the numerator of —k'ip) 

In order to factorize the numerator we introduce 

71(6; 6), 



-) 



72(6; 6) ^o(4-ff)' 

where kr = {kj^ — CiY^^ and 



(A.l) 



lo = hm — 7- = A11A33 ~ A 



13- 



(A.2) 



The function Ko{£,i ; ^2) thus defined tends to unity at infinity, is never zero and involves 
the ratio 71/72- The branch cuts of C2) and therefore ii'oCCi; C2) that fie in the lower 
half plane are presented in Fig. A.l. 



A Im ^1 



— Ki — le 
• ■ 

-K2 — ie 
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Re a 
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Fig. A.l. The branch cuts of function i<'o(Ci;C2) for a) ki imaginary and K2 real; b) ki and K2 real; c) 
Ki and K2 imaginary. 

Closing the integration contour in the lower half-space, we can write 
27ri t - a 



1 
27ri 



AlnA-oft;e2) 



dt. 



(A.3) 



where In Kq{£_i; ^2) is analytic in the upper half of the ^ plane, the jump over the cut is 
AlnAo(t;6)-lnA*(i,2li^)„h.A*(t,-2:i^^) (A.4) 

72(i;6) 72(i;6) 

and we have 

i^o(6; 6) = i^o+(ei; 6)ifo"(a; 6), (A.5) 



with 
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^o"(a;6) = <(-a;6) (A.6) 

analytic in the lower half of the ^i— plane. By changing t to —t, the formula (A. 3) can 
be simplified to 

27ri t + ^1 

When —^1 lies on the branch cut [k2, Ki], integration of the singular integrand can be 
avoided by using the formula 

6) = ^ ^ Ko{iv,i2). (A.8) 

^0 (,-?i;?2j 

Appendix B. Factorization of the denominator of Hi{4i;$2) 

In order to factorize the denominator of /ii(^i) we introduce the function 

V ■(: \ - -fc*((i Ill^iiM^ iziiiwii) . 

72(6;6) <i73(6;6) 

where we have 

a^oo 73(Ci;6) '- 4B1A33 J 

The function tends to unity at infinity, is never zero and involves the ratio 71/72- Using 
the same considerations as above, leads us to 

27ri t + ^1 

where the jump over the cut is 

Alni^i(t;6) = lnKr(t, ^^i^) - h,K{(t-^-^^^\ (B.4) 
We note that 

6) = 6)^r(a; 6), (b.5) 

with the function 

ifr(6;6) = i^i+(-a;6) (b.6) 

analytic in the lower half of the ^1— plane. It follows that when —^1 lies on the branch 
cut [k;2, integration of the singular integrand can be avoided by using the formula 

^t^i^-^i^) = ^/ -^ ^AirAi)- (B.7) 
^1 (-?i;?2j 
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Appendix C. Magnitudes of various diffraction coefficients for the axis of 
symmetry perpendicular to the crack 

In Fig. C.4 below, when the mode of the incident wave is qSV and the corresponding 
diffracted wave vectors are complex, we artificially set the diffraction coefficients to zero. 
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Fig. C.l. Magnitude D'^i^'IP'' versus 4, <p^'^'= = 90° and S'"" = Q°(a.), 30°(h), 60"(c), 90°(d), 1200(e), 
1500(f). 
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Fig. C.4. Magnitude D*'?^'''^^) versus = 60" and j^™'^ = 0°(a), 30"(b), 60''(c), 90'^(d), 1200(e), 

1500(f). 
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Appendix D. Magnitudes of various diffraction coefficients for the axis of 
symmetry lying in the crack plane perpendicularly to the crack edge 




Appendix E. The effect of cuspidal edges and conical points 



As mentioned at the end of Section 8, the cuspidal edges and conical points on the qSV 
wave surfaces lead to extra singularities in (123). In Fig. E.l we plot the "diffraction coef- 
ficients" as defined by (123) for the case of the symmetry axis perpendicular to the crack. 

In Fig. E.2 the corresponding quantities arc multiplied by ^27ri|sin d\S,^^\kf^^ ^ ; — fc™) 
to exclude the influence of the cuspidal edges. In Fig. E.3 further multiplication by 
A^"'(k'*'-''-'^) is carried out to exclude the influence of the conical points too. The remain- 
ing singularities are due to shadow boundaries. 
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Fig. E.l. Magnitude d(9^^'«'P) versus 4, <^'"= = 90° and i?'"'^ = 00(a), 30° (b), 600(c), 900(d), 1200(e), 
1500(f). 



41 




42 




43 



